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ABSTRACT 


An  analysis  of  steady  state  operation  of  a  vertically 
oriented,  variable  conductance  heat  pipe  is  presented.  The 
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I .  INTRODUCTION 


A.  BACKGROUND 

A  heat  pipe  is  a  closed,  evacuated  chamber  that  has  the 
inside  walls  lined  with  a  wick  structure,  saturates  /ith  a 
volatile  working  fluid.  The  steady  state  operation  inv  Ives 
vapor  heat  transfer  and  capillary  action.  A  comprehensive 
coverage  of  the  theory  and  operation  of  heat  pipes  can  be 
found  in  Dunn  and  Reay  [Ref.  11.  Briefly,  the  operation  of 
a  vertically  oriented  heat  pipe  will  be  discussed.  The  heat 
pipe  has  three  sections:  an  evaporator  section,  an  adiabatic 
section,  and  a  condenser  sc  ’lion.  The  heat  pipe  with  the 
three  sections  marked  is  shown  in  Figure  1.1.  The  lower  evapora 
tor  end  of  the  heat  pipe  is  heated,  causing  the  working  liquid 
in  the  wick  to  vaporize.  This  phase  change  causes  the  working 
fluid  to  absorb  latent  heat  and  the  vapor  pressure  increases 
in  the  evaporator  end.  The  working  fluid  is  drawn  away  from 
the  wick,  causing  the  liquid  pressure  to  decrease.  The  vapor 
moves  to  the  cooler  condenser  where  it  is  condensed  and  the 
working  fluid  gives  up  its  latent  heat.  The  condensate  is 
deposited  on  the  wick  which  causes  the  liquid  level  to  in¬ 
crease  thus  increasing  the  liquid  pressure.  Capillary  forces 
take  the  higher  pressure  liquid  from  the  condenser  to  the  lower 
pressure  liquid  end  of  the  evaporator.  The  capillary  forces 
are  assisted  by  gravitational  forces  when  the  heat  pipe  is 


vertically  oriented  as  shown  in  Figure  1.1.  It  is  important  to 
note  that  for  heat  pipes  the  evaporator  position  is  not  re¬ 
stricted,  but  since  the  experimental  data  and  numerical  model 
use  a  vertically  oriented  heat  pipe,  this  orientation  will  be 
used  throughout  this  thesis. 

The  great  advantage  of  heat  pipes  is  that  the  vaporization 
and  condensation  process  allow  for  high  heat  transfer  rates 
with  small  temperature  differences.  Two  other  advantages  are 
that  the  basic  heat  pipe  has  no  moving  parts  and  can  be  used 
in  many  geometrical  configurations.  The  heat  pipe  has  limi¬ 
tations  and  again,  a  comprehensive  presentation  of  these  limi¬ 
tations  can  be  found  in  Dunn  and  Reay  [Ref.  1].  Briefly,  the 
five  heat  pipe  limitations  are  as  follows: 

(a)  The  viscous  limit  is  basically  the  low  temperature  limit 
This  would  be  the  minimum  temperature  for  the  heat 

pipe  to  begin  operation.  Generally,  the  viscous  limit 
is  only  important  during  start  up  of  some  heat  pipes, 
i.e.,  liquid  metal  pipes. 

(b)  The  sonic  limit  describes  the  limit  of  the  vapor  flow 
velocity  at  the  heat  pipe  evaporator  exit.  The  re¬ 
sulting  choked  flow  condition  limits  the  axial  heat 
flux  capacity  of  the  heat  pipe. 

(c)  The  entrainment  limitation  is  due  to  the  removing  of 
the  liquid  before  it  can  return  to  the  evaporator  by 
vapor  flowing  at  high  velocity. 

(d)  The  wicking  or  capillary  limit  reflects  the  fact  that 
each  combination  of  wick  and  fluid  has  a  maximum 
capillary  head.  As  stated  earlier,  the  capillary  head 
must  increase  with  the  liquid  and  vapor  pressure  drops 
and  can  do  so  only  up  to  the  maximum  capillary  head. 

(e)  The  last  limit  is  the  boiling  limit.  If  temperatures 
in  the  evaporator  get  too  high  the  working  fluid  will 
begin  to  boil  within  the  wick.  This  will  disrupt  liquid 
flow  to  the  evaporator,  causing  the  heat  pipe  to  stop 
working . 


These  five  heat  pipe  limitations  vary  in  importance,  depending 
on  the  temperature  and  corresponding  axial  heat  flux  operating 
point.  Since  at  different  temperatures  one  can  encounter 
different  heat  pipe  controlling  limitations,  it  is  important 
to  have  a  basic  understanding  of  all  five  limitations. 

B.  VARIABLE  CONDUCTANCE  HEAT  PIPE 
1 .  Basic  Operation 

Variable  conductance  heat  pipes  are  described  in  detail 
in  Marcus  [Ref.  2] .  As  described  in  the  background  section,  a 
change  in  the  heat  load  will  automatically  change  the  operating 
temperature  and  pressure  of  the  heat  pipe,  given  a  constant  sink 
temperature.  This  is  an  undesirable  quality  for  many  situations 
where  it  is  desired  to  have  a  stable  source  temperature  under 
conditions  of  varying  heat  load  or  sink  temperatures.  An 
example  of  this  case  would  be  the  cooling  of  sensitive  elec¬ 
tronic  equipment.  Basically,  a  noncondensible  gas  is  introduced 
into  the  heat  pipe  along  with  the  working  fluid,  making  the 
heat  pipe  become  a  passively  controlled  device.  The  nonconden¬ 
sible  gas  occupies  a  certain  portion  of  the  condenser  section. 


depending  on  operating  conditions.  Increasing  operating  tempera 
ture  in  turn  increases  the  vapor  pressure  of  the  working  fluid 
which  compresses  the  noncondensible  gas  into  a  smaller  volume. 
This  provides  a  greater  active  condenser  area,  thus  allowing 
a  more  stable  source  temperature  than  a  conventional  heat  pipe. 
When  the  operating  temperature  decreases,  working  fluid  vapor 


pressure  falls  and  the  noncondensible  gas  expands,  blocking  a 
larger  portion  of  the  condenser.  The  result  is  a  passively 
controlled  variable  condenser  area,  which  in  turn  makes  a 
more  stable  source  temperature. 

2 .  Variable  Conductance  Heat  Pipe  Theory 

Marcus  [Ref.  2]  describes  a  flat  front  theory  and  a 
'diffuse  front  theory'  for  variable  conductance  heat  pipes. 

Flat  front  theory  assumes  the  noncondensible  gas  is  in  the 
form  of  a  plug  with  sharply  defined  vapor-gas  interface  which 
is  perpendicular  to  the  axis  of  the  pipe.  Since  this  model 
neglected  the  axial  conduction  in  the  heat  pipe  wall  and  diffu¬ 
sion  of  vapor  into  the  stagnant  gas  plug,  experimental  results 
[Ref.  2]  have  shown  this  model  to  be  inaccurate.  The  early 
model  of  flat  front  theory  provided  initial  insight  into  the 
operating  characteristics  of  gas  loaded  heat  pipes.  With  the 
incorporation  of  binary  diffusion  between  the  working  fluid 
vapor  and  the  noncondensible  gas  and  axial  conduction  in  the 
pipe  wall,  Marcus'  [Ref.  2]  'diffuse  front  theory'  establishes 
a  model  of  gas  loaded  heat  pipes  which  provides  better  simi¬ 
larity  between  computer  program  and  laboratory  results. 

Gravitational  effects  on  the  vapor-gas  interface  are  not 
included  in  flat  front  theory  or  the  'diffuse  front  theory. ' 
Gravitational  effects  can  greatly  distort  the  axial  temperature 
profile  when  the  working  fluid  and  noncondensible  gas  have 
significantly  different  molecular  weights.  This  fact  is  demon¬ 
strated  by  experimental  results  reported  by  Kelleher  [Ref.  3] . 


Batts  [Ref.  4]  conducted  variable  conductance  heat  pipe  experi 
ments  using  working  fluid/noncondensible  gas  combinations  of 
methanol/helium  and  methanol/krypton.  Kane  [Ref.  5]  developed 
an  analytical  model  of  variable  conductance  heat  pipes  which 
includes  the  effects  of  binary  mass  diffusion,  axial  pipe-wall 
conduction,  and  gravitational  effects  which  exist  at  the  vapor 
gas  interface  region  as  shown  by  Batts  [Ref.  4]  experiments. 


II.  FORMULATION  OF  AN  ANALYTICAL  MODEL 


A.  PROBLEM  STATEMENT 

The  purpose  of  this  chapter  is  to  formulate  the  equations 
describing  the  steady  state  operation  of  a  vertically  oriented, 
gas  loaded,  variable  conductance  heat  pipe.  As  in  Figure  1.1, 
the  condenser  is  above  the  evaporator  for  this  orientation. 

This  analysis  includes  the  effects  of  binary  mass  diffusion, 
axial  pipe  wall  conduction,  and  gravitational  effects  due  to 
the  difference  in  molecular  weight  between  the  noncondensible 
gas  and  the  working  fluid. 


B.  ASSUMPTIONS 

The  following  assumptions  are  used  to  develop  the  system 
of  equations: 

1.  Heat,  mass  and  momentum  transport  processes  are  treated 
as  one-dimensional  in  the  axial  direction  within  the 
vapor  space  of  a  vertically  oriented  heat  pipe. 

Kelleher  [Ref.  3]  shows  that  in  any  other  orientation, 
gravitational  effects  cause  highly  three-dimensional 
spatial  variations  when  the  molecular  weights  of  the 
working  fluid  and  noncondensible  gas  differ  substantially. 

2.  The  noncondensible  gas  acts  as  a  perfect  gas.  The 
mixture  of  noncondensible  gas  and  working  fluid  are 
treated  as  thermodynamically  independent  substances. 

3.  Because  of  direct  contact  between  the  vapor  and  the 
liquid  saturated  wick,  the  working  fluid  vapor  is  at 
saturation  conditions  at  all  times. 

4.  The  thermal  resistance  of  the  wick-wall  structure,  in 
the  radial  direction,  is  negligible  in  comparison  with 
the  thermal  resistance  due  to  convection  outside  the 
pipe. 


C.  ANALYSIS 


Kane  [Ref.  5]  presents  a  proposed  solution  for  this  prob¬ 
lem.  The  following  is  a  significant  change  to  his  proposed 
solution  with  modifications  being  noted  as  they  occur  in 
this  section.  Variables  with  the  subscript  'a'  apply  to  the 
working  fluid  and  variables  with  the  subscript  'b'  apply  to 
the  noncondensible  gas. 

Utilizing  the  assumption  that  the  noncondensible  gas  is 
a  perfect  gas,  the  equation  for  noncondensible  gas  pressure 
is ; 


P 


b 


p  —  T 


(II.  1) 


The  principles  of  mixtures  of  thermodynamically  independent 
substances,  which  are  found  in  Reynolds  and  Perkins  [Ref.  6], 
are  used  to  define  the  total  density  and  the  total  pressure. 


P 


Pa  ^ 


(II. 2) 


P  =  P  +  P. 

a  b 


(II. 3) 


The  definition  of  mass  fraction  is  used  to  form  the  equation 


for  the  working  fluid  vapor  mass  fraction. 


\ 


The  solution,  utilizing  the  conservation  of  mass,  is 
based  on  a  differential  element  representing  a  section  of 
the  heat  pipe  condenser.  Figure  2.1.  The  conservation  of 


I 

\ 

I 

I 

1 


mass  of  species 


'a'  for  the  differential  element  yields; 


‘•'a'^a  S'0a''a>'^’‘'*c  "aV* 


This  equation  sets  the  mass  of  species  'a'  into  the  control 
volume  equal  to  the  mass  of  species  *a'  out,  including  a 
term  to  account  for  the  removal  of  species  'a'  by  condensation 
The  process  of  condensation  is  accounted  for  by  multiplying 
the  condensation  mass  flux  times  the  area,  defined  by  the 
vapor  space  perimeter  times  the  length,  dx,  of  the  differen¬ 
tial  element.  After  basic  algebra,  the  equation  for  conser¬ 
vation  of  mass  of  species  'a'  becomes; 


0  = 


m  P 
a  V 


^  "-c  dS^^'^a^a^ 


By  using  Pick's  law  of  binary  diffusion  [Ref. 7],  and  solving 
for  p  V  ,  the  resulting  formula  is; 

di  3 


P  V 
^a  a 


dw 


The  final  form  for  the  conservation  of  mass  for  species  'a' 
becomes ; 

.  P  ,  ,  dw 

"'a  =  0  . 


(II. 5) 


since  species  'b'  is  a  noncondensible  gas,  the  condensation 
mass  flux  of  species  'b'  equals  zero.  The  resulting  equation 
for  the  conservation  of  mass  for  species  'b'  is: 


b 


pV) 


A 

- 


0  . 


The  mass  continuity  of  the  mixture  is  derived  by  adding 
the  conservation  of  mass  of  species  'a'  with  the  conservation 
of  mass  of  species  'b' . 


*  ^  H  A 

"‘a  ■  d5r^P°(-dr  “dST 

c 


)  = 


d  d 

Because  w  +ui,  is  equal  to  one,  and  3— (w  +aj,  )  =  xT(i)  =  0, 
a  D  ^  dx  a  t)  dx 

the  final  result  of  the  mass  continuity  of  the  mixture  is: 


c 


(II. 6) 


Figure  2.2  is  a  differential  element  of  the  vapor  space  in 
the  condenser  and  is  used  to  derive  the  conservation  of  linear 
momentum  equation  for  the  mixture  of  noncondensible  gas  and 
working  fluid  vapor.  The  conservation  of  linear  momentum  for 
the  mixture  of  noncondensible  gas  and  working  fluid  vapor  is 
derived  by  setting  the  summation  of  external  forces  acting  on 
the  control  volume  in  the  x  direction  equal  to  the  net  efflux 
of  linear  momentum. 


I  =  /  V(pV  dA) 


The  summation  of  external  forces  acting  in  the  x  direction  are 


p  gA^  dx  -  T  P  dx  -  dp  A 
c  V  c 


These  terms  represent  the  force  due  to  gravity  acting  on  the 
control  volume  minus  the  forces  due  to  wall  shear  stress  and 
pressure  acting  on  the  control  volume.  The  wall  shear  stress 
can  be  calculated  by: 


T 


V 


2 


where  A  is  a  dimensionless  friction  factor.  If  developed 
laminar  flow  is  assumed  the  Darcy  friction  factor,  A  = 
could  be  used  with  the  Reynolds  number  defined  as 


P  VD 


Re  = 


The  equation  for  wall  shear  stress  then  becomes: 


T 


where  y  is  the  viscosity  of  methanol  vapor,  V  is  the  velocity 
of  the  vapor  and  is  the  diameter  of  the  vapor  space.  The 
terms  to  represent  the  net  efflux  of  linear  momentum  are: 


/  V ( pvdA) 

CS  ,x 


■  2 

V  m  P  dx  +  (pv  +d (pv) (v+dv) ) A  -  pV  A 

3  cl  V  L*  C 

V  m  P  dx  +  A  pvdv  +  A  vd(pv)  . 

^  3.  V  C  ^  C 


By  setting  the  suimnation  of  external  forces  in  the  x  direction 
equal  to  the  net  efflux  of  linear  momentum  and  combining  like 
terms,  the  result  is: 


dP 

dx 


d  ,  2, 

pg  -  33^(pv  ) 


8pv 

D 

V 


(II. 7) 


The  above  equation  represents  conservation  of  linear  momentum 
for  the  mixture  of  noncondensible  gas  and  working  fluid  vapor. 
Kane  [Ref.  5]  assumes  zero  velocity  for  the  noncondensible  gas 
and  derives  an  equation  for  conservation  of  linear  momentum 
of  species  'b. '  The  species  'b*  conservation  of  linear  momen¬ 
tum  equation  is  not  necessary  to  derive  the  model  in  this 
analysis.  Since  the  conservation  of  linear  momentum  for  the 
mixture  of  the  noncondensible  gas  and  working  fluid  vapor  is 
the  more  consistent  approach.  Equation  (II. 7)  is  used  through¬ 
out  this  thesis. 

Figure  2.3  is  a  differential  element  of  the  condenser  wall 
and  represents  the  wick  and  pipe  wall  combination  as  one  ele¬ 
ment.  The  wick-wall  combination  is  treated  like  a  fin  and  is 
used  to  derive  the  equation  for  the  conservation  of  energy. 
Heat  transfer  is  represented  by  the  latent  heat  deposited  on 
the  'fin'  by  condensation  of  the  working  fluid  vapor  plus  the 
heat  transfer  from  Fourier's  law  of  conduction  entering  the 


Differential  Element  for  Energy  Balance 


Heat  transfer  out  is  represented  by  Fourier's  law  of  conduction 
leaving  the  'fin'  plus  the  energy  which  is  convected  away  from 
the  fin. 


Q 


out 


Q  +  dQ  +  h  P  (T  -  T  )dx 
C  C  S  “ 


The  equation  for  conservation  of  energy  is  derived  by  setting 
heat  transfer  in  equal  to  heat  transfer  out  of  the  'fin.' 


d^T 

dx2 


hP 

< 


■(T  -T 


h.  P  . 
kA  a 


=  0 


(II. 8) 


The  experimental  data  from  Batts  [Ref.  4]  is  used  to 
correlate  this  analysis.  The  working  fluid  is  methanol, 

CH^OH.  The  working  fluid  vapor  in  the  heat  pipe  is  in  direct 
contact  with  the  working  fluid  liquid  in  the  wick  and  is 
therefore  considered  to  be  a  saturated  vapor  at  the  local 
temperature  of  the  wick-wall  combination.  Eubank  [Ref.  8] 
tabulates  data  for  the  saturation  properties  of  methanol. 

These  tables  are  used  to  generate  functions  of  local  temperature 
for  partial  pressure,  partial  density  and  heat  of  vaporization 
of  saturated  methanol  vapor. 

(a)  The  partial  pressure  of  saturated  methanol  vapor  is  a 
function.  Fa,  of  local  temperature. 


P,  =  Fa(T) 


(II. 9) 


(b)  The  partial  density  of  saturated  working  fluid  vapor 


III.  NUMERICAL  MODEL  FORMULATION 

A.  VARIABLE  ELIMINATION 

The  system  of  eleven  equations  can  ultimately  be  reduced 
to  a  system  of  three  equations  by  eliminating  variables. 

Equation  (II. 9)  is  eliminated  by  using  the  equation  of  state, 
Fa(T),  for  P^.  The  partial  density  of  the  saturated  working 
fluid  vapor,  (Equation  II. 10)  is  eliminated  by  using  Sa(T). 
Hjg(T),  Equation  (II. 11),  is  used  to  eliminate  The  non¬ 

condensible  gas  pressure,  Pj^,  is  eliminated  by  using  the  perfect 
gas  equation  of  state  (Equation  II. 1).  The  treatment  of  worki  ig 
fluid  vapor  and  noncondensible  gas  as  thermodynamically  inde¬ 
pendent  substances  gives  total  pressure  as  the  sum  of  partial 
pressures  and  total  density  as  the  stun  of  partial  densities. 
Total  pressure  is  eliminated  by  using  the  sum  of  partial  pres¬ 
sures,  Equation  (II. 3).  The  noncondensible  gas  density  is 
eliminated  by  solving  Equation  (II. 2)  for  Pj^  and  replacing  Pj^ 
with  p  -p  .  The  definition  for  mass  fraction  of  the  working 
fluid  vapor.  Equation  (II.4),  is  used  to  replace  total  density. 

After  making  the  preceding  seven  substitutions,  four  equa¬ 
tions  with  four  unknowns  are  left.  The  four  remaining  depen¬ 
dent  variables  which  are  functions  of  axial  position  are: 

■  2 

(a)  m  ,  working  fluid  condensation  flux  [Ibm/ft  sec] 

3 

(b)  V,  velocity  of  the  mixture  [ft/sec] 

(c)  (jj  ,  working  fluid  mass  fraction  [0  <  w  <  1] 

(d)  T,  temperature  [°R]. 


The  four  remaining  equations  are: 


c  a 


=  0 


(III.l) 


"a  JT  *  35'=^  V) 
c 


^(SaD 

dx  '  u)  dx  ' 
a 


=  0 


(III. 2) 


P  Tfi  P 

(§3:  g,  .  d  (§a  „2  _  _!!.  v_i  ._v  8>i.  V 
“a  "a  *o  “a  °v 


(III. 3) 


(d) 


dfr 

dx^ 


P  P  . 

— ^  h(T  -T  )  +  ^  m 

KA  KA  ™a 

P  P 


(III. 4) 


Furthermore,  Equation  (III.l)  can  be  eliminated  by  solving  this 
equation  for  working  fluid  condensation  flux,  m  ,  and  substi- 

a 

tuting  that  into  the  three  remaining  equations. 

The  resulting  three  equations  with  dependent  variables  of 
velocity,  working  fluid  mass  fraction,  and  temperature  are: 


(a)  -^(—  V)  +  ^(Sa  V)  -  ^(— D-3-^)  =  0 

dx  (i)  dx  dx  0)  dx 

a  a 


(III. 5) 


cl  D 


ax  o)  0)  dx  oj 


Pv8y 
A  D 


V 


(III. 6) 


(C) 


KA 


h  (T  -T  ) 


“fq  C  ^,Sa 

KA  dx'o) 

P  a 


V)  = 


(III. 7) 


This  set  of  three  equations  and  three  unknowns  describe  the 
operation  of  a  vertically  oriented, gas  loaded, variable  conduc¬ 
tance  heat  pipe. 

B.  BOUNDARY  CONDITIONS 

These  equations  describe  the  thermodynamics  of  the  heat 
pipe  condenser.  The  independent  variable,  the  distance  along 
the  heat  pipe,  is  measured  with  the  origin  at  the  condenser 
end  farthest  from  the  evaporator.  The  total  condenser  length 
is  defined  as  L.  Based  on  the  system  of  three  equations,  two 
boundary  conditions  are  required  for  temperature  and  mass  frac 
tion  and  one  boundary  condition  is  required  for  velocity. 

From  experimental  data,  the  temperature  at  the  far  end  of 
the  condenser  and  in  the  evaporator  are  known.  The  boundary 
conditions  for  the  dependent  variable  of  temperature  are; 

when  X  =  0,  T  =  T  ,  where  T  equals  the  temperature  at  the 

far  end^of  the  condenser. 

when  X  =  L,  T  =  Tg,  where  T^  equals  the  temperature  at 

the  end^of  the  condenser  closest  to 
the  evaporator. 

The  highest  velocity  of  the  working  fluid  vapor  is  reached 


in  the  evaporator  exit.  This  velocity  can  be  calculated  from 
the  known  heat  input  to  the  heat  pipe.  The  vapor  velocity  in 
the  evaporator,  V  ,  is  equal  to  the  heat  transfer  rate  into 


I.%«l  ia*  !«»  »■>>  #■* ‘fi*'!^*  li*  tA*  IW*  tail  Wf 1^’ia*  WlW'ta*  W  M  tk* 
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the  evaporator  divided  by  the  product  of  local  density  of  the 
working  fluid  vapor,  cross  sectional  area  of  the  vapor  space, 
and  the  local  heat  of  vaporization  of  the  working  fluid.  Due 
to  the  physical  boundary  at  the  end  of  the  pipe,  where  x  =  0, 
the  velocity  at  the  farthest  end  of  the  condenser  away  from  the 
evaporator  is  equal  to  zero.  The  boundary  conditions  for  the 
dependent  variable  of  velocity  are: 


when  X  =  0 ,  V  =  0 . 0  and 


when  X  =  L,  V  =  V 


P  A  h^ 
^a^  c  fg^ 
e  ^e 


Based  on  the  system  of  three  equations,  only  one  of  the  above 
boundary  conditions  is  required.  The  boundary  condition  at  the 
origin  is  used  for  the  iteration  scheme  and  is  used  to 
normalize  velocity. 

The  temperature  and  pressure  at  the  farthest  end  of  the 
condenser  from  the  evaporator  are  known  from  laboratory  data. 
Assuming  thermodynamically  independent  substances,  the  total 
pressure  is  the  sum  of  partial  pressures  and  the  total  density 
is  the  sum  of  partial  densities.  The  saturation  pressure  of 
the  working  fluid  vapor,  P, ,  could  be  determined  from  tabulated 
data  by  entering  the  saturation  pressure  table  with  condenser 
temperature.  An  example  of  this  table  can  be  found  in  Eubank 
[Ref.  8]  for  methanol  saturated  vapor.  The  noncondensible  gas 
pressure,  Pj^,  can  be  calculated  by  subtracting  the  working 
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fluid  vapor  pressure  from  the  total  pressure.  With  this  data,  j 

the  initial  mass  fractions  can  be  established  by  using  the  | 

following  procedure  from  Reynolds  and  Perkins  [Ref.  6];  ‘  j 

1 


total 


total 


X 


M 


a  a 


Lbm  _ 

Lbmole  mixture 


for  species  'a' 


Xv^M.  =  — s - ! — ■  -  _  for  species  'b' 

^b  b  Lbmole  mixture 

^a^a 

"=0  '  >‘a”a  ^XfaMb 

'=0  ■  Xa^a  ■"Xb'V 

At  the  end  of  the  condenser  closest  to  the  evaporator,  it  is 
assumed  that  the  concentration  of  noncondensible  gas  is  zero 
so  that  at  this  point  u)  is  equal  to  one.  The  boundary  condi- 
tions  for  the  dependent  variable  of  mass  fraction  for  the 
working  fluid  vapor  are; 

when  X  =  0 ,  o)  =  ui  ,  and 
a  Sq' 


when  X 


L. 


0) 


1.0  . 


C.  EQUATION  NONDIMENSIONALIZATION 


To  generalize  the  equations  and  facilitate  their  solution 
it  is  convenient  to  use  dimensionless  variables  defined  as 
follows . 

(a)  The  independent  variable,  condenser  length,  is  defined 
by; 

X 

=  L 

where  x  is  the  distance  along  the  heat  pipe,  measured  from 
the  condenser  end  farthest  from  the  evaporator,  and  L  is  the 
total  condenser  length. 

(b)  A  dimensionless  temperature  can  be  defined  by: 


T  -T 

00 


where  T  is  the  local  temperature,  T^  is  the  evaporator  tempera 
ture,  and  T^  is  the  ambient  temperature.  In  this  analysis,  th 
first  and  second  derivative  of  temperature  with  respect  to  dis 
tance  are  required. 

dT  _  ^  ^  ^  _  "^e  '^<”d6 

dx  ~  d0  dn  dx  ~  L  dn 

and 


dx'dx^ 


T  -T 

d  .  e  °°  d9  >  dn 
dn '  L  dn  dx 


(c)  The  second  dependent  variable,  working  fluid  mass 
fraction,  is  dimensionless  by  definition  and  its  numerical 
value  is  greater  than  or  equal  to  zero  and  less  than  or  equal 
to  one.  The  first  derivative  of  working  fluid  mass  fraction 
with  respect  to  distance  is: 

doj,  doj  j  ,  do) 

_ a  _  a  ^  _  1  a 

dn  dn  dx  L  dn 

(d)  Velocity  of  the  working  fluid  vapor  is  made  dimension¬ 
less  by  normalizing  local  vapor  velocity  with  the  evaporator 
vapor  velocity,  V^.  The  nondimensional  velocity  and  the  first 
derivative  of  velocity  with  respect  to  distance  are  then 
defined  by: 


V* 


V 

V 

e 


dV  ^  dV  dV*  ^  ^  ^  dV* 

dx  dv*  dn  dx  L  dn  ’ 


(e)  In  order  to  complete  nondimensionalizing  process,  the 
following  other  dimensionless  variables  are  used. 

(i)  The  working  fluid  vapor  density  is  normalized  with 


respect  to  the  evaporator  vapor  density. 


(ii)  The  working  fluid  vapor  pressure  is  compared 


to  the  evaporator  vapor  pressure. 


Fa* 


(iii)  The  working  fluid  heat  of  vaporization  is  normalized 
with  respect  to  the  condenser  working  fluid  vapor  heat  of 
condensation. 


H 

H 


f2_ 

fg,c 


(iv)  The  working  fluid  vapor  viscosity  is  normalized 
with  respect  to  the  evaporator  vapor  viscosity. 


y* 


After  substitution  of  nondimensional  terms,  the  three 
equations  which  describe  operation  of  a  vertical  variable 
conductance  heat  pipe  become; 


(a) 


d  .Sa*-V* 

dn '  u) 

ci 


f^(Sa.V) 


D  ,d  Sa*  *^^a. 
LV^Mn  dn 


0  (III. 8) 


The  dimensionless  coefficient  from  Equation  (III. 8)  is 


e  e 


Lg.  Sa*. 
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b  e 


d  ,Sa*V*‘',  .  „*  d  ,Sa*V*, 

arr^-n; — ^  ■"  ^ 

cl  ci 


P  8vi  L 
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A  D  Sa  V 
c  V  e  e 
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(III. 9) 


The  dimensionless  coefficients  from  Equation  (III. 9)  are 


N2  = 


Sa  V 
e  e 
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N3  =  - -  ,  N4  = 


b  e 


T  -T  ' 
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P  hL 
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V  e 

K  D„Sa^ 
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The  two  nondimensional  terms  of  interest  from  Equation  (III.  10) 


P  hL 

N7  =  r —  and  N8 

I\  A 
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A  Sa„  L 
fg,c  c  e  e 

K  A  (T  -T^) 
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Numerical  evaluation  of  the  dimensionless  coefficients  is 
required  to  determine  the  relative  significance  of  each  term. 


Two  experimental  runs  from  Batts  [Ref.  4]  are  used  to  evalu¬ 
ate  the  dimensionless  coefficients.  A  low  powered  run 
(Q  =  22.9  watts)  with  methanol  as  the  working  fluid,  and  krypton 
as  the  noncondensible  gas,  is  chosen.  This  run  demonstrates 
an  example  of  the  heavier  molecular  weight  noncondensible  gas 
diffusing  through  a  large  portion  of  the  condenser.  A  second 
higher  powered  run  (Q  =  70  watts)  with  methanol  as  the  working 
fluid,  and  helium  as  the  noncondensible  gas,  is  chosen  because 
it  demonstrates  an  example  of  the  lighter  molecular  weight 
noncondensible  gas  occupying  a  small  portion  of  the  condenser. 
The  convection  heat  transfer  coefficient  for  the  heat  pipe  is 

^5  1 

set  equal  to  1.0  - x — ,  representing  free  convection.  A 

hr.ft.  °R 

system  thermal  conductivity  of  approximately  4.56  - ~ — 

hr.ft.  °R 

results  from  a  parallel  circuit  combination  of  conductive 
thermal  resistances  for  the  stainless  steel  pipe  wall  and 
methanol  filled  stainless  steel  wire  mesh  wick.  Kane  [Ref.  5] 

calculated  representative  binary  mass  diffusion  coefficients 

2  2 
of  0.4  ft  /hr  for  krypton-methanol  and  1.9  ft  /hr  for  helium- 

methanol.  These  values  for  binary  mass  diffusion  are  also 

used  in  this  analysis.  The  experimental  values  that  are  used 

to  calculate  the  dimensionless  constants  are  listed  in  Table 

I.  The  values  of  the  dimensionless  constants  for  the  two 

experimental  runs  are  listed  in  Table  II. 

Based  on  the  tabulated  numerical  results.  Equation  (III. 9) 


is  modified  to: 


TABLE  I 


NUMERICAL 

Constants  for  both  runs : 

R  =  1545.33  ft  Lbf/Lbmole  °R 

g  =  32.2  ft/sec^ 

L  =  3.5  ft 

A  =  .0199  ft^ 
c 

A  =  .001917  ft^ 

P 

D  =  .15917  ft 

V 

P  =  .524  ft 
s 

=  .500  ft 

h  =  1.0  Btu/hr  ft^ “R 

K  =  4 . 5585  Btu/hr  ft  °R 

Low  powered  krypton/methanol 

Q  =  22.9  watts  =  78.1348  Btu/hr 

T  =  592°R 
e 

T  =  543°R 
c 

T  =  536°R 

OO 

D  =  0.4  ft^/hr 

M.  =  8  3.8  Lbm/Lbmole 


DATA 

Gas  constant 

Acceleration  of  gravity 

Length  of  condenser 

Cross  sectional  area  of 
vapor  space 

Axial  cross  sectional  area 
of  wick-wall  combination 

Diameter  of  vapor  space 

Periphery  of  outside  of 
pipe 

Periphery  of  vapor  space 

Convection  heat  transfer 
coefficient 

Thermal  conductivity 


TABLE  I  (CONTINUED) 


=  .041339  ft/sec 
Fa^  =  10.24  Lbf/in^ 

=  500.21  Btu/Lbm 

fg,c 

Sa^  =  0.054871  Lbm/ft^ 

=  224.502  xio"^  Lbf  sec/ft^ 

High  powered  helium/methanol 

Q  =  70  watts  =  238.84  Btu/hr 

T  =  640°R 
e 

T  =  540°R 
c 

T  =  536°R 
00 

D  =  1.9  ft^/hr 

Mj^  =  4.003  Lbm/Lbinole 

V  =  .050485  ft/sec 
e 

Fa^  =  28.52  Lbf/in^ 

=  501.11  Btu/Lbm 

fg,c 

Sa^  =  0.14446  Lbm/ft^ 
u  =  243.993  xlO"^  Lbf  sec/ft^ 


TABLE  II 


DIMENSIONLESS  CONSTANTS 


Low  Powered  Krypton/Methanol  High  Powered  Helium/Methanol 


LV 


^  =  7.679  X  lo"'^ 
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Sa  V 
e  e 
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M. 
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T  -T 
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V  e 


A  D  Sa  V 
c  V  e  e 
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P  hL 
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K  A 


=  734.553 
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N2  = 


=  19.4581  X 10  N3  = 


N4  = 


N6  = 


N7  = 


N8  = 


Pa 


Sa 

e  e 


=  359.166  X lo' 


R(T  -T  ) 
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M. 

b  e 


=  507.224  X 10 


e  ° 


5.154 


N5  = 

V 
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.=  44218.0 


P  8  y  L 
V  e 


A  D  Sa  V 
c  V  e  e 


=  4.76143 


P  hL 
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K  A 


=  734.553 


A  Sa  V  L 
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K  Ap(T^-T^) 


581.362 


1008.30 
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e  e 


M.  V 
b  e 
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The  viscous  term  and  the  two  inertia  terms  are  neglected  be- 

1  3  6 

cause  their  order  of  magnitude  is  10  compared  to  10  and  10 
for  the  remaining  terms.  The  nondimensional  equations  (III. 8), 
(III. 10),  and  (III. 11)  are  used  to  develop  the  analytical  model 
of  the  vertically  oriented, variable  conductance  heat  pipe. 


D.  FORMULATION  OF  THE  NUMERICAL  MODEL 

Equation  (III.  8)  is  integrated  and  is  solved  for  the  first 
derivative  of  the  working  fluid  mass  fraction.  Integrating 
both  sides  of  Equation  (III.  8)  gives: 


Sa*V* 


+  Sa’^'V*  -  Nl 


Sa’ 


duj 


U) 


dn 


=  c 


(III. 12) 


The  lower  endpoint  is  used  (n  =  0)  to  evaluate  the  constant 
of  integration.  When  n  =  0,  V*  at  the  end  of  the  pipe  equals 
zero.  The  constant  of  integration  then  becomes: 
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-Nl 


dw 
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(III .13) 


By  replacing  the  value  for  the  constant  of  integration  and 
solving  for  doj  /dn ,  the  final  version  of  Equation  (III.  12)  is: 


N1  Nl  Sa*  Sa^  dn  q 


(III. 14) 


The  derivative  in  Equation  (III. 11)  is  expanded  to  get  the 


following. 


+N4)) 


(e  .N4)  +  sa*  ^  ^(9  +N4) 
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The  chain  rule  is  applied  to  . 


dSa* 

dn 


dSa*  dSa  ^  ^ 
dSa  dT  d0  dn 


T  -T 
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dFs* 

Applying  the  same  procedure  to  ,  the  result  is: 


dFa*  _  /^e  dFa  d0 

dn  ~  '  Fa  '  dT  dn 
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(III. 17) 


By  combining  terms  and  solving  for  the  final  form  for 

equation  (III. 11)  is: 


(a)  The  coefficient  of  ^  is  set  equal  to  G(a)^,0) 


G  = 


dFa  .  »3(T^-T,)  1 


Fa  dT  Sa  dT  'w 


(—  -  1)  (0  +N4) 


+  N3  (Sa*)  (—  -1)  . 
“a 


(111.18) 


d  0 

(b)  Solving  for  the  result  is 


d0  rNSSa*  .  N3Sa* 


OJ 


(III. 19) 


Equation  (III.  10)  is  integrated  and  the  result  is: 


d0 

dn 


-N7 


0dr| 


m  H|^ 


Sa*V* 


u> 


C  . 


(III. 20) 


Using  the  boundary  condition  of  n  =  0  where  V*  =  0.0,  the 
constant  of  integration  is: 


C 


(III. 21) 


The  final  version  of  Equation  (III.  10)  which  is  solved  for  V* 
is : 


V*  =  ( 


d0 

arr 


-N7  (  f  0dn )  w,  -  0), 

cL  '  a  a 


dn 


)/(N8H*gSa*) 


(III. 22) 
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The  final  system  of  three  equations  with  nondimensional 
variables  is: 


(a)  Equation  (III. 14),  -^(V* ,a)^ ,-^|  ); 


-V*  0,^  P(,  du.^ 

N1  N1  Sa^  dn  j 


jQ  du) 

(b)  Equation  (III. 19),  g^(w^,0,-^); 


de  _  rN5Sa*  ,  N3Sa*  a  ...a,  ,  / 

dTT  -dlf(e+N4)]/G(u)^,e) 

a  0) 

a 


where 


N2(T  -T  )  N3(T  -T  )  ,„  , 

/-/..  AX  _  e  «'  dFa  .  '  e  »'  dSa ,  1  ix  x 

G(w^,9)  -  ^ - ^ - g_(_ -1)  (0+N4) 

e  e  a  ■ 


+  N3(Sa*)  (—  -1) 
“a 


(c)  Equation  (III. 22),  V*  (w  ,  /'^0dn,^,^  ) 

a  0  an  an  Q 


The  boundary  conditions  for  this  set  of  equations  are  the 


following: 


(a)  For  nondimensional  temperature: 


n  =  0.0 


0 


C 


n  =  1.0 


0  =  1.0 


The  constant  a  is  not  known  but,  based  on  laboratory  data,  can 
be  estimated  and  through  an  iteration  process  the  correct  value 
can  be  obtained.  A  proposed  iteration  scheme  is  presented  in 
Chapter  IV. 

(b)  For  working  fluid  vapor  mass  fraction; 


n  =  0.0 

a  a^ 

n  =  1.0  0),  =  1.0 

d 

doj 

The  constant  0  is  not  known,  but  can  also  be  estimated  for  an 
initial  value  and,  through  an  iteration  scheme,  a  correct  value 
can  be  obtained. 

(c)  For  working  fluid  vapor  dimensionless  velocity: 


E.  AUXILIARY  EQUATIONS 

A  means  to  derive  saturated  working  fluid  heat  of  vapori¬ 
zation,  density,  and  pressure  is  required  before  a  numerical 
solution  is  attempted.  Using  methanol  as  the  working  fluid  and 
data  from  Eubank  [Ref.  8]  for  saturation  conditions,  the 
following  functions  are  derived. 

1 .  Polynomial  for  the  Working  Fluid  Heat  of  Vaporization 
Tabulated  values  of  the  heat  of  vaporization  for  methanol 

are  plotted  on  Figure  3.1  with  a  solid  line.  A  Fortran  program, 
polyfit,  is  used  to  obtain  the  coefficients  of  the  following 
polynomial : 

(T)  =  aT^  +  bT  +  c  (III. 23) 

fg 

where  a  =  -9.96  xio”^  Btu/Lbm ( ®R) ^ ,  b  =  .730  Btu/Lbm(“R),  and 
c  =  397.99  Btu/Lbm.  The  values  of  H^g(T)  using  this  polynomial 
are  plotted  on  Figure  3.1  with  a  dotted  line. 

2 .  Polynomial  for  the  Working  Fluid  Vapor  Saturation 
Density 

Tabulated  values  for  methanol  saturated  density  are 
plotted  on  Figure  3.2  with  a  solid  line.  Again  the  Fortran 
program,  polyfit,  is  used  with  the  following  modification; 

T*  =  -  where  T  ^  =  500°R  and  ^ 

T  -  ret 

ref 


A  polynomial  in  the  variable  C  is  used  to  describe  the 
saturation  density. 


The  polynomial,  Sa(C).  is  then; 


Sa{C)  =  aC^  +  +  cC  +  d 


(III. 24) 


where  a  =  .1293  Lbm/ft^,  b  =  -1.0475  Lbm/ft^,  and 
c  =  2.8868  Lbm/ft^  and  d  =  -2.7037  Lbm/ft^.  The  value  of 
dSa/dT  is  also  required  for  the  numerical  solution.  This  is: 


^  2bC  + 


(III. 25) 


3 .  Polynomial  for  the  Working  Fluid  Vapor  Saturation 
Pressure 

Tabulated  values  for  methanol  saturation  pressure  are 
plotted  on  Figure  3.3  with  a  solid  line.  Polyfit  is  used  with 
the  same  modif ication . 

The  polynomial  Fa(C)  is: 


Fa(C)  =  +  cC  +  d 


where  a  =  24.5994  Lbf/in^ ,  b  =  -191.2472  Lbf/in^, 
c  =  498.4940  Lbf/in^,  and  d  =  -439.9883  Lbf/in^.  The  first 
derivative  of  Fa(C)  with  respect  to  temperature  is  equal 


^  =  (3a +  2bC^  +  cC)/T_„^ 


(III. 26) 


IV.  DISCUSSION  AND  RECOMMENDATIONS 


A.  DISCUSSION 

The  final  version  of  the  model  for  a  variable  conductance 
heat  pipe  with  a  vertical  orientation  is  a  system  of  three 
ordinary,  non-linear,  coupled,  differential  equations.  These 
three  equations  are  (III. 14),  (111.19),  and  (III. 22). 

The  remaining  challenge  for  the  numerical  model  is  to  pro¬ 
pose  an  iteration  scheme  that  can  solve  this  system  of  three  equa¬ 
tions.  A  special  thank  you  goes  to  the  collaboration  of  Dr.  M.D. 
Kelleher  and  Dr.  K.T.  Yang  for  the  following  iteration  technique. 
The  steps  of  the  iteration  procedure  are  as  follows: 

(a)  To  start  the  iteration  process,  assume  a  form  for  the 
relationship  between  0  and  such  that  o)^  =  Cj^0  +  C2 .  This 
linear  relationship  is  reasonable  because: 

(i)  When  ri  =  0 ,  both  u  and  0  are  values  close  to 

zero. 

(ii)  When  ri  =  L  both  w  and  0  are  equal  to  1.0. 

(iii)  This  step  is  required  for  the  first  iteration 

only. 

The  equations  for  this  step  are  then: 


“a  ~  *“1®  ^2 


(IV. 1) 


(IV. 2) 
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This  establishes 


Then,  using  Equation  (III. 19),  solve  for 
a  first  degree  equation,  in  terms  of  0  only  which  can  then 

be  integrated  numerically  to  establish  a  first  iteration 
relationship  for  9(ri). 

dfi 

(b)  With  this  relationship  for  6  and  ^  defined  from  step 
(a)  above.  Equation  (III. 22)  can  be  solved  for  V*. 

(c)  Once  step  (b)  is  completed.  Equation  (III.  14)  can  be 

solved  for  dui  /dn .  The  integral  of  dw  /dn  is  then  used  to 
3i  3 

define  (jj  .  The  process  can  then  be  repeated. 

In  summary,  the  first  iteration  goes  through  steps  a-c 
and  the  remaining  iterations  go  through  steps  b-c  until  there 
are  no  further  changes  in  the  form  of  0 ,  oj  ,  and  V*.  In 

3i 

theory,  the  outcome  of  this  numerical  iteration  is  a  table 
which  gives  dimensionless  temperature,  0;  velocity,  V*;  and 
mass  fraction  of  the  working  fluid  vapor,  co  ,  versus  non- 
dimensional  condenser  length,  n . 


B .  RECOMMENDAT IONS 

1.  This  numerical  analysis  was  written  to  yield  results  via 
a  computer  program.  The  Continuous  System  Modeling 
Program  III  [Ref.  9]  is  a  computer  package  that  could 

be  used  to  solve  the  system  of  equations.  The  CSMP 
[Ref.  9]  includes  a  powerful  integral  statement  which 
carries  out  a  Fourth  Order  Runge-Kutta  Method. 

2.  This  model  assumes  constant  values  for  the  binary  mass 
diffusion  coefficient  and  the  convection  heat  transfer 
coefficient.  Once  numerical  values  are  obtained  from 
the  iteration  scheme,  the  accuracy  can  be  improved  by 
using  analytical  expressions  instead  of  constants  for 
the  binary  mass  diffusion  coefficient  [Ref.  10]  and  the 
convection  heat  transfer  coefficient  [Ref.  11]. 

3.  An  attempt  could  be  made  to  model  the  variable  conduc¬ 
tance  heat  pipe  that  is  not  in  a  vertical  orientation. 
Kelleher  [Ref.  3]  shows  that  a  stratification  of  the 


working  fluid  and  noncondensible  gas  occurs  when  the 
heat  pipe  is  in  the  horizontal  position.  This  position 
makes  the  transport  processes  inside  the  pipe  three 
dimensional,  which  would  require  a  finite  difference 
computer  model  for  solution. 
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